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ABSTRACT 

We report new upper limits on any possible long-term time variation of the 
ratio of the electron to proton masses, ( m e /m p ), the fine-structure constant, (a), 
and the quantity a 2 g p (m e /m p ), where g p is the proton gyromagnetic ratio. These 
limits are based on extremely high precision observations of H 2 , Si 3+ , C° and H° 
in high-redshift quasar absorption lines. They amount to 95% confidence ranges 
of (—7.6 —> 9.7) x 10 -14 yr _1 for ( m e /m p ), (—4.6 —■> 4.2) x 10“ 14 yr _1 for a and 
(—2.2 —> 4.2) x 10- 15 yr _1 for a 2 g p (m e /m p ), where the elapsed time has been 
computed for a cosmology with H 0 = 75 km s^ 1 Mpc^ 1 and qo = 0.5. 


Subject headings: atomic data — cosmology:observations — cosmology:theory — 
quasars:absorption lines 


1. Introduction 


Current laboratory, geophysical and astrophysical constraints exclude any significant 
evolution of most of the dimensionless physical constants over cosmological time (Sisterna & 
Vucetich 1990). However, refining these limits remains important as they provide constraints 
on theories such as Kaluza-Klcin and superstrings which allow solutions in which there can 
be time variation of these constants. We report improved astrophysical constraints on the 
fine-structure constant, a and the ratio of the electron to proton masses, (m e /m p ). 

Savedoff (1956) first pointed out the possibility of using differential measurements of 
redshifted lines in cosmologically distant objects to test the evolution of certain of the di¬ 
mensionless physical constants; since then various measurements have been outlined and the 
observational constraints greatly improved. Recent summaries of these astrophysical mea¬ 
surements may be found in Potekhin & Varshalovich (1994) and Levshakov (1994). In this 
paper we describe the use of extremely high signal-to-noise and high resolution (R — 36, 000 
to 47,000) spectra of quasar absorption lines obtained with the Keck 10 m telescope to 
tighten the evolutionary constraints on a, ( m e /m v ) and (a 2 g p m e /m p ) by roughly an order 
of magnitude over previous measurements. Note that, in translating limits to time varia¬ 
tion rates we have assumed a flat cosmology (qo = 0.5, A = 0) with a Hubble constant of 
H 0 = 75 km s -1 Mpc -1 . 


2. Constraints on ( m e /m p ) 


The presence of cosmological evolution in ( m e jm v ) can be tested by using observations 
of molecular hydrogen in quasar absorption line systems, assuming that all other dimen¬ 
sional constants, particularly the fine-structure constant, do not vary in such a way as to 
compensate for the variation in (■ m e /m p ). This was first pointed out by Thompson (1975), 
who noted the different dependence of the electronic, vibrational and rotational energy levels 
on (■ m e /m p ). There is, however, a general unavailability of data at high redshift: to date, 
only one quasar, PKS 0528 — 250, is believed to have molecular hydrogen in its spectrum, at 
a redshift of 2.811 (Chaffee, Foltz & Black 1988). Varshalovich and Levshakov (1993) have 
analyzed this data to obtain alu limit of d(m e j m p )/dt < 3 x 10~ 13 yr _1 , assuming a time 
difference of 1.3 x 10 10 yr between z = 2.811 and z = 0. 

We have re-observed PKS 0528 — 250 at much higher spectral resolution (R = 36, 000) 
and signal-to-noise, obtaining a 5.4 hour spectrum with the Keck Telescope which unam¬ 
biguously confirms Chaffee et al.’s identification of molecular hydrogen, as demonstrated in 
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Fig. 1 which shows the 0-0 and 4-0 bands of the Lyman series. An analysis of the physical 
properties of the molecular hydrogen itself is given elsewhere (Songaila & Cowie 1995). 

Table 1 summarises the properties of those lines that are isolated molecular hydrogen 
lines, based on their measured column densities and FWHM. In the Born-Oppenheimer 
approximation, the molecular hydrogen energy levels can be written, 

E = E el + 'j E mb + (^) E rot . (1) 

\rripj \ m pJ 

So, the energy shift in any vibration-rotation transition j in the Lyman series has the form, 

AB, = a d + b, (A) + c, (A) (2) 

and the difference in energy between two transitions is, 

(A Ej ~ AEi) ~ bji (^) + ° ji (^) ‘ (3) 
To lowest order, a change in ( m e /m p ) = fx induces a change in A Ei — A Ej such that 

5 ± ~ ( 2AE 1 ) (4) 

H ~ c \AEi-AEj) 1 ; 

where Sv is the mean offset compared to the laboratory value, of the energy difference 
between the two sets of lines, when that offset is represented as a velocity difference. While 
the above description is in fact adequate for the analysis and gives a clear insight into the 
method, we have in practice used the exact relationship, 

7 - »• - « (?) - ». (?) • 

where the k coefficients have been tabulated by Varshalovich & Levshakov (1993). 

A regression analysis based on the data of Table 1 then gives a best fit value Sfi/n — 
8 x 1CA 5 and a 95% confidence range of —5.5 to 7 x 10~ 4 . The r.rn.s. velocity spread is 
2.5 km s _1 . For our ‘standard’ cosmology, the time to z — 2.811 is 7.2 x 10 9 yr and our 
95% confidence limit for jx/fx is —7.6 to 9.7 x 10 -14 yr -1 . This is approximately an order of 
magnitude improvement over the previous result. 
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3. Constraints on a 


The separation between the wavelength (Ai) corresponding to the transition 2 S \/2 —> 2 
P 3/2 and that (A 2 ) corresponding to the transition 2 S ' 1 / 2 —» 2 P 1/2 in an alkaline ion is pro¬ 
portional to a 2 to lowest order in a (Bethe & Salpeter 1977). Writing A = (2/3)Ai + A 2 , 


(Ai — A 2 ) 
A 


AA 2 
—r - ~ a. 


( 6 ) 


So any change in a will result in a corresponding change in AA in the separation of the 
doublets in a high- 2 ; quasar (as was first used by Bahcall & Salpeter [1965] and Wolfe, 
Brown & Roberts [1976]) 


5a 

a 



S( AA) 
2AA 



( 7 ) 

( 8 ) 


The most recent analyses using this method (Potekhin & Varshalovich 1994) have used very 
large inhomogeneous samples of C IV, N V, O VI, Mg II, AI III and Si IV doublets and 
obtained a 95% confidence limit of \a~ 1 (da/dz)\ < 5.6 x 10“ 4 . 


We have improved this analysis using a homogeneous sample of doublets observed at 
high spectral reolution. However, before proceeding further, it is important to note that, 
for a given accuracy in determining A v/c, the sensitivity of determining A a is inversely 
proportional to the relative splitting of the doublet (AA/A). For this reason, Si IV, the most 
widely spaced of the doublets (AA1393.755, 1402.770 A; Toresson 1960) is by far the most 
sensitive probe of a, followed by AI III and Mg II. We have therefore restricted our analysis 
to the Si IV doublet. Choosing a doublet with large (AA)/A also alleviates a second problem, 
namely that, for many of these doublets, the laboratory wavelength separation is not known 
to very high precision. This enters directly in the ( 8v/c ) term as a systematic uncertainty 
which is minimised for large (AA)/A. For Si IV, Martin and Zalubas (1983) estimate the 
uncertainty in Toresson’s (1960) wavelengths at a maximum of 5 mA, which, if adopted 
as an extreme upper bound to the uncertainty in the doublet spacing, translates into an 
uncertainty of 1.0 km s _1 in An, and currently imposes a fundamental limit of 2.7 x 10~ 4 on 
the accuracy to which 5a/a can be determined. 


We have identified a number of Si IV doublets in high- 2 ; quasar absorption-line systems 
observed with the Keck HIRES spectrograph at R = 36, 000. The redshift systems used in 
the analysis were seen not only in Si IV but also in lines of other species, such as C IV and 
Ly a; therefore, the identifications do not depend solely on the spacing of the Si IV lines. 
Table 2 lists the relative velocity shifts for Si IV. As is illustrated in Fig. 2, the offsets have 
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been measured by cross-correlating the two members of the doublet. The r.m.s. uncertainty 
in this procedure is 1.1 km s _1 , determined from measurements of the more common C IV 
doublet. The 95% confidence limits inferred from Table 2 are ( Sa/a ) = (—2.2,1.6) x 10~ 4 , 
comparable to the systematic error from the uncertainty in the laboratory wavelengths. 
Combining the two errors in quadrature, we have \5a/a\ < 3.5 x 10 -4 and our results imply 
\oT l da/dz\ < 1.1 x 10 -4 , a factor of 5 lower than Potekhin and Varshalovich’s (1994) result. 
Further improvements will require either new laboratory determinations of the Si IV doublet 
separation, or local observations of the Si IV doublet spacing in the nearby interstellar and 
intergalactic gas. The later can be carried out by the spectrographs on the Hubble Space 
Telescope. To be useful for this purpose, high resolution measurements with the Goddard 
High Resolution Spectrograph would be required. Current published values are not adequate. 


4. Constraints on a 2 g p (m e /m p ) 


The ratio of the frequencies of the hyperhne 21 cm absorption transition of neutral 
hydrogen (u a ) to an optical resonance transition (z/&) has the dependence 


Va 


Vb 


ot g p 


m P 


m ri 


= X 


(9) 


so that evolution of this quantity will result in a difference in the measured redshifts of 21 cm 
and optical absorption 

fSx\ 

Sz = z opt - z 2 i = (1 + z) ( — j . (10) 

The current best limits on this quantity are given by Tubbs and Wolfe (1980), who found 
(1 /x)(dx/dz) < 1.1 x 10” 4 for absorption in the quasar Q1331+170. The redshift of the 
hyperhne absorption in Q1331+170 is known to very high precision (z 2 i = 1.77642±2 x 10 -5 ; 
Wolfe & Davis 1979) and the uncertainty in the Tubbs and Wolfe (1980) estimate arises 
primarily in the imprecision of the optical redshift. Recently, Songaila et al. (1994) have 
presented high S/N observations of C° absorption and fine-structure C°, finding z opt = 
1.77644 ± 2 x 10“ 5 . Because the C° arises in the same cloud components responsible for the 
21 cm absorption, the comparison of z 2 1 and z opt should be relatively secure against different 
kinematic structures being present in the two measurements. Combining the two redshifts, 
we find 

A Y 

( 11 ) 


— = 7x 10“ 6 ±1.1 x KT 5 
X 

corresponding to a 95% confidence range for (1 /x)(dx/dz) of (—2.2,4.2) x 10~ 15 yr~ 4 for our 
standard cosmology. 
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5. Conclusion 


The three constraints presented here can be compared with existing laboratory and 
astrophysical constraints. The limits on a 2 g p (m e /m p ) currently constitute the tightest as- 
trophysical constraint on the evolution of these dimensionless parameters, implying \a/a\ < 
2.1 x 10 -15 yr -1 and \(m e /m p ) / (m e /rn p )\ < 4.2 x 10 -15 yr -1 at the 95% confidence level, in 
the absence of relative cancellation. For a/a, the astrophysical measurement is now compa¬ 
rable to the best local tests (Sisterna & Vucetich 1990), which give \a/a\ < 1.3 x 10 -15 yr^ 1 
(95% confidence), but, because of the larger timeline, may provide more stringent constraints 
on models with non-linear evolution of a. 


We would like to thank Art Wolfe for many very helpful discussions. The authors 
were visiting astronomers at the W. M. Keck Observatory, which is jointly operated by the 
California Institute of Technology and the University of California. 
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Table 1. Molecular Hydrogen Lines in PKS 0528 — 250 


Line 

X a 

/ 'vac 

k h 

(< 8v/c ) (xlO 5 ) c 

0-0 R(0) 

1108.128 

7.79 

0.1 

0-0 R(l) 

1108.633 

8.25 

-0.9 

1-0 R(0) 

1092.195 

0.72 

-0.6 

1-0 R(l) 

1092.732 

1.22 

0.2 

1-0 P91) 

1094.052 

2.38 

0.7 

1-0 R(2) 

1094.225 

2.62 

1.2 

1-0 P(2) 

1096.438 

4.54 

1.3 

2-0 R(l) 

1077.699 

-5.24 

-1.1 

2-0 P(l) 

1078.925 

-4.18 

0.3 

2-0 R(2) 

1079.226 

-3.81 

0.2 

2-0 P(2) 

1081.226 

-2.05 

-0.9 

3-0 P(l) 

1064.605 

-10.21 

0.5 

3-0 R(2) 

1064.995 

-9.73 

-0.4 

3-0 P(2) 

1066.900 

-8.10 

-1.2 

4-0 R(0) 

1049.367 

-17.27 

0.8 

4-0 R(l) 

1049.959 

-16.67 

-1.0 

4-0 P(l) 

1051.032 

-15.76 

-0.5 

4-0 R(2) 

1051.498 

-15.18 

1.4 

4-0 P(2) 

1053.284 

-13.67 

0.0 


a Dabrowski 1984 

hPotekliin & Varshalovich 1994 

c Frame of reference is such that ( 8v/c) = 0 
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Table 2. Si IV Doublet Separations 


Quasar 

z 

Sv 

0302-003 

2.785 

0.8 

0528-250 

2.813 

0.5 

0528-250 

2.810 

0.4 

0528-250 

2.672 

-2.1 

1206+119 

3.021 

-0.5 

2000-330 

3.551 

-1.5 

2000-330 

3.548 

1.1 

2000-330 

3.332 

2.3 

2000-330 

3.191 

-2.2 
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Fig. 1.— Spectra of the 0-0 and 4-0 vibration transitions of the Lyman series of molecular 
hydrogen at z — 2.811 in the quasar PKS 0528 — 250. The spectrum is shown in the rest 
frame with the dashed lines indicating the wavelengths of the J < 2 transitions in each band. 
The remaining (generally broader) lines in the spectra are primarily Lyman alpha forest lines 
of intergalactic neutral hydrogen. 


Fig. 2.— Cross correlation of the two lines of the Si IV doublet for a cloud in the spectrum 
of the quasar Q0302 — 003. The quantity (1 — S'(A)) is plotted versus v, where S is the 
normalised spectrum for the stronger line, and the quantity 2(1 — S( A)) (dashed line) is 
overplotted versus v — Av for the weaker line. Here Av is the measured offset. The vertical 
dashed lines show the velocity range over which the cross correlation function was computed. 












